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Abstract 

We calculate the semileptonic form factors f+^ n (q 2 ) and f+~* v (q 2 ) from QCD sum rules 
on the light-cone (LCSRs), to NLO in QCD, and for small to moderate q 2 , < q 2 < 
16 GeV 2 . We include in particular the so-called singlet contribution, i.e. weak annihilation 
of the B meson with the emission of two gluons which, thanks to the U(1)a anomaly, 
couple directly to rf 1 '. This effect is included to leading-twist accuracy. This contribution 
has been neglected in previous calculations of the form factors from LCSRs. We find that 
the singlet contribution to f+~* v can be up to 20%, while that to f+~^ v is, as expected, 
much smaller and below 3%. We also suggest to measure the ratio B(B — > v(ev)jB{B — > 
rjeu) to better constrain the size of the singlet contribution. 
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Figure 1: Flavour-singlet contribution to a generic B — * rf transition. 

1 Introduction 

B — > rf'' transitions are interesting for a number of reasons: at tree-level, they involve a 
b —>■ u transition and hence are sensitive to the CKM matrix element \V u b\. Its precise de- 
termination is crucial for the interpretation of the "tension" [1] that has emerged between 
the determination of \ V u b\ from, on the one hand, inclusive semileptonic B — > X u lv decays 
[2], and, on the other hand, global fits [1, 3] and the exclusive decay B — > niv [4, 5, 6, 7]. 
The inclusive value of \V u b\ is larger than that from other determinations and hints at a 
non-zero new-physics contribution to the B d mixing phase 0^, i.e. 4>d ^ 2(3 [8]. While an 
analysis of all available experimental and theoretical information on B — > i\iv found no 
"significant" disagreement between the exclusive and the inclusive values of \V u b\ [6], the 
situation has changed very recently, when the HPQCD lattice collaboration reported a 
mistake in their calculation of the form factor f+^ n published in Ref. [7]; the corrected 
form factor is larger and hence yields a smaller \V u b\ [9]. The authors of Ref. [6] have 
since then published an update [10] of their previous analysis and now conclude that the 
exclusive value of \V u b\ is in perfect agreement with the determination from global fits 
and that "the hints of a disagreement with inclusive determinations of \V u b\ are strength- 
ened". Also very recently, Neubert has argued [11] that the value of \V u b\ obtained by 
the HFAG collaboration [12] is dominated by observables with small efficiency and that, 
selecting observables with maximum efficiency instead, the resulting \V u b\ is smaller than 
the HFAG average. Given this situation it is important to collect information on \V u b\ 
also from other exclusive processes. B — > rf'^lv decays offer the opportunity for doing so. 

Another reason why B — ► rf 1 ' transitions are interesting is their sensitivity to r\-rf 
mixing and the effects of the U(1)a anomaly, which is responsible for the large mass 
of the rf and also induces potentially large flavour-singlet contributions to amplitudes 
involving rf''. Indeed the unexpectedly large branching fractions of inclusive B — > r/'X and 
exclusive B — > rfK decays, as compared to e.g. B — ► tc transitions, have been attributed 
to an enhanced flavour- singlet contribution [13], which is defined as the amplitude for 
producing either a quark-antiquark pair in a singlet state (uu + dd + ss) which does not 
contain the 5's spectator quark, or a pair of gluons, followed by hadronization into an 
rf'\ A generic contribution of this type is shown in Fig. 1. In Ref. [14] it was found that 
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a rather large singlet-contribution of ca. 30% to the form factor would bring the 

central values of theoretical predictions for B — > rj'K observables in QCD factorisation 
into good agreement with experimental results, although the theoretical uncertainties are 
too large to allow a definite conclusion on the size of the singlet contributions. On the 
other hand, a more recent analysis of B decays with isosinglet final states, formulated 
in SCET, finds that, because of large experimental uncertainties of the data used to fit 
non-perturbative parameters, the singlet contribution to form factors is consistent with 

[15]- 

While the interplay of singlet and octet contributions is well understood at the level 
of local matrix elements, i.e. decay constants (wave functions at the origin) [16, 17, 18], 
less is known about the shape of these wave functions, which are relevant for dynamical 
quantities like form factors. In frameworks based on QCD factorisation the mesons' 
Fock-state wave functions enter in the form of light-cone distribution amplitudes (DAs). 
Constraints on the leading parameters of these DAs have been obtained from the analysis 
of the 77^7 transition form factor [19, 20, 21] and of the inclusive decay Y(1S) — > rj'X [21]. 
In principle, these DAs can also be constrained from a measurement of the form factors 
of B — > rf'\ for instance from B(B — > rj'£u)/B(B — > r/lis), as suggested in Ref. [22]. 

Despite the strong phenomenological interest in the size of the singlet contribution 
to f+~* v , there is, to the best of our knowledge, only a single calculation available, 
based on the perturbative QCD approach [23]. Ref. [23] finds that this contribution is 
negligible in f+^ v , and reaches a few percent in f+^ v . Another well-known method 
for the calculation of B — > light meson form factors are QCD sum rules on the light 
cone (LCSRs) [24, 25, 26]. Ref. [26], for instance, provides form factors for B — > {n,K,rj) 
decays, but does not include the singlet contribution to B — > 77, nor a calculation of B — > 77' 
form factors. It is the purpose of this paper to remedy this situation and complete the 
calculation of B — > light pseudoscalar meson form factors from LCSRs by including also 
the flavour- singlet contributions. 

Our paper is organised as follows: in Sec. 2 we define the two most common 77-7/ 
mixing schemes and review rf'^ DAs. In Sec. 3 we derive LCSRs for the B — > rf'^ form 
factors. In Sec. 4 we present results and conclude. 

2 77 and 77' Mixing and Distribution Amplitudes 

There are two different mixing schemes in use to describe the 77-77' system: the singlet- 
octet (SO) and the quark-flavour scheme (QF) [16]. In the former, the couplings of the 
relevant axial-vector currents to the meson P = 77, 7/ are given by 

(0\JU\ p (p))-ifpP, (< = 1,8), (1) 

where J^ 5 denotes the SU(3) F -octet and JK the SU(3) F -singlet axial-vector current, re- 
spectively. The four parameters f p define two decay constants /j of a hypothetical pure 
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singlet or octet state \r)i) and also two mixing angles 6>j via 

f" fi \ - ( cos6 * ff 8 o 



ft ft) \sm9 8 cos^AO/i '' (2) 

The advantage of this scheme is that the impact of the U(1)a anomaly is plainly localised 
in f ± , via the divergence of the singlet current J^ 5 , while Q { ^ and / 8 ^ f w are SU(3)p- 
breaking effects. By the same token, the SO scheme also diagonalises the renormalisation- 
scale dependence of parameters and hence is very useful for checking the cancellation 
of divergences in perturbative calculations: / 8 and Q; L are scale-independent, while fi 
renormalises multiplicatively [27]: 

In the QF mixing scheme, on the other hand, the basic axial-vector currents are 

J% = "y= {ui^U + d^ 5 d) , J* 5 = S7m75S , (4) 

and the corresponding couplings to P = rj, rj' are given by 

(0\JU\P(p)) = ifpP, (r = q,s). (5) 
In complete correspondence to (2) one has 




cos0 9 -sin0 s \ / f g 
sin (f) q cos 4> s ) V f s 



(6) 



The basic difference to the SO scheme is that now the difference between the two angles (f) qjS 
is not caused by SU(3)f effects, like that between 6\ and 6*8, but by an OZI-rule violating 
contribution, as explained in Ref. [17]. While the numerical values of 9i differ largely, with 
typical values 6 S pa —20° and 6 X pa —5°, one finds 4> s — <p q <, 5°, with <fi q pa <fi s pa 40° [16, 17]. 
This led the authors of Ref. [16] to suggest the QF scheme as an approximation to describe 
rj-rj' mixing, based on neglecting the difference (p q — (p s (and all other OZI-breaking effects): 

= <frq,s, (p q -(l) s =0. (7) 

The advantage of this scheme is that it has only 3 parameters, f q , f s and 0, which implies 
that the mixing of states is the same as that of the decay constants: 

1] \ _ f cos (f) - sin <f) \ f r] q \ ^ 



rj J \ sin <p cos <p j \r\ 



Is 



The disadvantage is that, due to the neglection of OZI-breaking effects, the renormalisa- 
tion-scale dependence of f\ is not reproduced - as it is induced precisely by OZI-breaking 
terms [17]. While this is not really an issue numerically, as the scale-dependence of f\ is 
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a two- loop effect, Eq. (3), the problem of the incompatibility of the QF scheme with the 
scale-dependence of parameters will come back at the level of non-local matrix elements, 
i.e. DAs, see below. 

Given enough data to fix all independent parameters, there is no reason to prefer the 
QF over the SO scheme. For DAs, however, the SO scheme leads to a proliferation of 
unknown parameters, while the QF scheme is more restrictive, see below. For this reason 
we decide to use the QF scheme in this paper. Its basic parameters have been determined 
as [16] 

/, = (1.07 ± 0.02)/,, f s = (1.34 ± 0.06)/, , = 39.3° ± 1.0° . (9) 
This can be translated into values for the SO parameters as 



/> \/^ + ^ = (l-26±0.04)/ w , 



h = \J\P q +\f 2 s= (1-17 ±0.03)/,, 

9 8 = <j)- axctan[V2 f s /f q } = -21.2° ± 1.6° , 

6 1 = <f>- fxctfm[V2 f q /f a ] = -9.2° ± 1.7° . (10) 

Note that in the QF scheme / 9)S are scale-independent parameters, and so is fi as obtained 
from the above relations. The SO decay constants can be expressed in terms of the QF 
ones and the angle as 



f 8 f 1 
fl 1 /?/ 




(11) 



Let us now turn to light-cone DAs, that is the extension of matrix elements like (1) 
and (5) to those over non-local operators on the light-cone. This paper is not the place 
to give a thorough discussion of the properties of DAs, for which we refer to reviews [28] 
and to Refs. [29, 30]. Suffice it to say that the DAs are ordered in terms of increasing 
twist, with the minimum, or leading, twist for meson DAs being two. Motivated by the 
structure of the evolution of DAs under a change of the renormalisation scale //, they are 
expanded in terms of so-called asymptotic DAs multiplied by Gegenbauer polynomials. 
In the context of this paper it is important to recall that the U(1)a anomaly induces, in 
addition to two-quark DAs, also two-gluon DAs, of both leading and higher twist. Some 
properties of these higher-twist DAs have been studied in Ref. [21]. In this paper we only 
include the effects of the leading-twist two-gluon DA, which is justified as its effects turn 
out to be small and higher-twist DAs are estimated to have even smaller impact. We will 
come back to that in Sec. 4. 
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We define the twist-2 two-quark DAs of rf'^ as [20] 

(0\^z)Crf'y fi [z,-z]*(-z)\P(p))=i(pz)r P fdue^^u). (12) 

Jo 

Here is a light-like vector, z 2 = 0, and [x,y] stands for the path-ordered gauge factor 
along the straight line connecting the points x and y, 



[x, y] = P exp 



igl dt ( x - y )^(tx + (l-t)y) 
. Jo 



(13) 



u (1 — u) is the momentum fraction carried by the quark (antiquark) in the meson, £ is 
short for 2u — 1. 02-p( M ) * s the twist-2 DA of the meson P with respect to the current 
whose flavour content is given by Cj, with ^ = (w, d, s) the triplet of light-quark fields in 
flavour space. For the SO currents, one has C\ — 1/ a/3 and C 8 = As/a/2, while for the QF 
currents C g = [\p!C\ +C%)/ a/3 and C s = (Ci — V2C 8 ) / y/3, with Aj the standard Gell-Mann 
matrices. 

The gluonic twist-2 DA is defined as 1 

(0\G, z (z)[z,-z]G^(-z)\P(p)) = \ (pzf^fh ^due^'m^u) . (14) 

In order to perform the calculation of the correlation function defined in the next section, 
we also need the matrix element of the meson P over two gluon fields. Dropping the 
gauge factor [z, —z], one has 

Because of the positive G-parity of r\ and ?/, the two-quark DAs are symmetric under 
u <-> 1 — u: 

<&. P {u) = 4p(l - u) ; (16) 
they are expanded in terms of Gegenbauer polynomials as 



J 2:P\ 



«)=6u(l-«) (l+ E ^(^^(Oj (< = l,8,(z, a ); (17) 

V n=2,4,... / 



the quark Gegenbauer moments. As for the two-gluon DAs, the asymptotic DA is 
it 2 - 7_1 (l — li) 2 - 7-1 with j = 3/2 the lowest conformal spin of the operator the expansion 
goes in terms of Gegenbauer polynomials Cn 2 ■ One can show that i\)\. p is antisymmetric: 

^p(u) = -Vf :P (i - «) ; (18) 



1 This definition refers to the "cr-rescaled" DA 0g in Ref. [20] with a = \/3/Cp. It agrees with that 
used in Refs. [21, 23], which means that we can use their results for the two-gluon Gegenbauer moment 
B\ without rescaling. 
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in particular du^-pi 11 ) = an d the local twist-2 matrix element {Q\G^ Z G^ Z \P) van 
ishes. The non-vanishing coupling (Q\G a pG a P\P) induced by the U(1)a anomaly is < 
twist-4 effect. The corresponding matrix elements are given, in the QF scheme, by [16]: 



(0\a s GG/(An)\r] q ) = f s {m 2 v - mj,) sin cos , 
(0\a s GG/(4ir)\r) s ) = j q {m% - mj)/v^sin0cos( 



(19) 



We will estimate the size of these effects in Sec. 4. There are no twist-3 two-gluon DAs and 
the remaining twist-4 DAs also have vanising normalisation, see Ref. [21]. The conformal 
expansion of the twist-2 two-gluon DA reads 



(20) 



n=2,4,. 



with the gluonic Gegenbauer moments B^ 9 . In this paper, we truncate both (f) l 2 . p and ip^-p 
at n = 2. This is due to the fact that our knowledge about these higher-order Gegenbauer 
moments is very restricted. An estimate of the effect of higher Gegenbauer moments in 
2;7r on the B — > n form factor has been given in Ref. [31], based on a certain class 
of models for the full DA beyond conformal expansion. The effect of neglecting o n >4 we 
found to be very small, ~ 2%. We expect the truncation error from neglecing B g n>A to be 
of similar size. 

(j)\. P and ip2-p m i x upon evolution in /i, see for instance Ref. [20]. This amounts to a 
mixing of a^' 1 and B^' 9 , resulting in the renormalisation-group equation, to LO accuracy, 



d 

^ dji 



An 



( 100 10 \ 

81 

- 36 22 



/ 



(21) 



where for simplicity we have dropped the superscript P. We only quote the solution for 



a\{n) 



49 \ 



+ 



2 2 V / 276T 
5 



J \2 



/(2A)) 



a\(no) 



9 V / 276T 



L 7 2 "/(2/3 ) _ £7+7(2/30) 



b 9 M 



(22) 



with L = a s (fj,) / a s (fj,o) and the anomalous dimensions 7^ = (149 ± v / 276l)/9. This is to 
be compared to the evolution of the octet Gegenbauer moment: 



(23) 



Numerically, the evolution of a\ does not differ much from that of a|, for a wide range 
of B 9 ,: assume a|(l GeV) = a\{l GeV), as is the case for a strict imposition of the QF 
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Figure 2: Dependence of 4(2.4 GeV) on Bf(lGeV), Eq. (22), for o^(lGeV) = 0.2. 

scheme. Choose af(lGeV) = 0.2, as indicated by our knowledge of twist- 2 DAs of the 
re; then we have 0^(2.4 GeV) = 0.137 from (23); 2.4 GeV is a typical scale in the calcu- 
lation of form factors from LCSRs. In Fig. 2 we show the results of the evolution of the 
singlet Gegenbauer moment a\ from 1 to 2.4 GeV, from Eq. (22), for the range of gluon 
Gegenbauer moments \B 2 (1 GeV)| < 10. Evidently the impact of the different anomalous 
dimensions of a\ and a 2 is negligible (a<j,(2.4 GeV) = 0.137 for B 2 = 0) and the mixing of 
B 2 into a\ is smaller than 20% within the range of B 2 considered. 

At this point we would like to come back to the impact of evolution on the consistency 
of the QF scheme. We introduce the twist-2 two-quark DAs 2 , % — 1, 8, q, s, corresponding 
to the basis states \r]i) in the SO and QF scheme, respectively. We then have, in terms of 
the quark valence Fock states \qq) and \ss) [20]: 

k> ~ 4>l(u)\qq) + <PT\u)\ss) , \ Vs ) ~ <j>° Z \u)\qq) + <j> s 2 (u)\ss) , (24) 

where qq is shorthand for (uu + dd) / \/2 and 

0« = I (0* + 20 2 ) , 0* = I(20 8 + 2 O, 0^ = ^(01-0*). (25) 

In the QF scheme, the "wrong-flavour" DA 0° ZI , which is generated by OZI-violating 
interactions, is set to 0. Once this is done at a certain scale, however, the different 
evolution of a\ and a^, Eqs. (22) and (23), will generate a non-zero 0!? ZI already to LO 
accuracy. A consistent implementation of the QF scheme hence requires one to either set 
a^' 8 = and also B-^ = 0, or to set a 8 = a\ and neglect the different scale-dependence 
of these parameters. In practice, however, the QF scheme is an approximation anyway, 
motivated by the observed smallness of one parameter, the difference of mixing angles 
S — <fi q . The induced non-zero DA 0° ZI is numerically very small for the scales relevant 
for our calculation, /j, — 1 GeV and 2.4 GeV. We hence implement the QF scheme for DAs 
as follows: we set (p\ = s , at the scale /i — 1 GeV, which, by virtue of (25), implies 02 = 0| 
at the same scale. We then evolve 02 according to the scaling-law for the octet Gegenbauer 
moment, Eq. (23). 2 We also set ^f;?? = V'f-rj'j a g a i n an y SU(3)f breaking of this relation is 

2 This is equivalent to imposing the QF-scheme relation a\ = a| as the scale [i — 2.4 GeV and defining 
B\ as B| (2.4 GeV). 
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expected to have only very small impact on ^ ( . The twist- 2 parameters used in our 
calculation are then reduced to 2: 0,2 and -Bf • For error estimates, we will also sometimes 
distinguish between a\ and a\ ■ 

As far as numerics is concerned, we assume that the bulk of SU(3)F-breaking effects 
is described by the decay constants via f q 7^ f n , and that SU(3)f breaking in Gegenbauer 
moments is subleading. This motivates setting a\ = , with a£(l GeV) = 0.25±0.15 as an 
average over a large number of calculations and fits to experimental data [30]; this number 
also agrees with a recent lattice determination [32] . justified discussed in 

Ref. [30], there is no evidence for noticeable SU(3)-breaking effects between a\ and af and 
the main SU(3)-breaking in the DAs is due to non-zero odd Gegenbauer moments. In this 
work we only need a|, and as a QCD sum rule for this parameter would look essentially the 
same as that for aj, except for a slightly different value for the decay constant, f n 7^ f q , 
and different numerical values for the continuum threshold s and the window in the Borel 
parameter M 2 , we see no plausible source for large SU(3) breaking between a\ and a\. 
To the best of our knowledge, no calculation of B 9 2 is available. Results from fits to data 
have been obtained from the 77^7 transition form factor, yielding B\ (1 GeV) = 9 ± 12 
[20], and the combined analysis of this form factor and the inclusive decay Y(1S) — > rj'X 
yielding 51(1.4 GeV) = 4.6 ± 2.5 [21]. These results, however, have to be taken cum 
grano salis as they are highly correlated with the simultaneous determination of a\ and 
a% from the same data, yielding a\(l GeV) = -0.08 ± 0.04, a^(lGeV) = -0.04 ± 0.04 
[20] and 4(1.4 GeV) = a\{lA GeV) = -0.054 ± 0.029 [21]. The same analysis applied to 
the 7T7 form factor returns a£(l GeV) = —0.06 ± 0.03 [33]. These results are not really 
compatible with those from the direct calculation of a\ from lattice and QCD sum rules; 
in particular the sign of a\ is unambiguously fixed as being positive. A possible reason for 
this discrepancy is the neglect of higher-order terms in the light-cone expansion and that, 
in addition, as one of the photons in the process is nearly real with virtuality q 2 ~ 0, one 
also has to take into account long-distance photon interactions, of order 1/ y^g 2 " [34]. For 
this reason, we assume the very conservative range -Bf (2.4 GeV) = 0±20 in the remainder 
of this paper. 

As far as higher-twist DAs are concerned, we only need those involving currents with 
flavour content qq = (uu + dd) / \/2. In line with the implementation of the QF scheme 
for twist- 2 DAs, we include SU(3)f breaking only via the decay constants and set 

l(0m z )C q [z,-z]r*(-z)\ri('\p)) = ^(0\d(z)[z,-z}Tu(-z)\n-(p)), 

J v (0 Jtt 

-^(0\^(z)[z,vz}G(vz)C q T[vz,-z]^(-z)\r ] ^(p)} = 
J v (>) 

(0\d(z)[z,vz]G(vz)F[vz, -z]u{-z)\ir-{p)) , (26) 

where T is the relevant Dirac structure and G(vz) the gluon field-strength tensor. The 
precise definitions of all twist-3 and 4 DAs, as well as up-to-date numerical values of 
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the 7r's hadronic parameters can be found in Ref. [30]. Let us shortly comment on the 
validity of this treatment for twist-3 two-quark DAs. As is well known, the normalisation 
of these DAs is given, for the ir, by f^m^/ {2m q ) and enters the light-cone sum rules for 
B — > 7r transitions as a l/m& correction, see explicit formulas for the corresponding D 
form factor in Ref. [35] . Although suppressed by one power of the heavy quark mass, this 
contribution is numerically non-negligible due to the chiral enhancement factor. Following 
the above implementation of SU(3) breaking, we set / 7r m^/(2m 9 ) — > f q m1/(2m q ) for r\ q 
(the corresponding quantity for r\ s is not needed). In contrast, the inclusion of all SU(3) 
effects leads one to consider the quantity 

K = f q (m 2 v cos 2 + m 2 , sin 2 0) - v^/^m 2 , - m 2 ) sin cos ; (27) 

the normalisation of the twist-3 DAs of i] q is given by h q /(2m q ). To leading order in 
the chiral expansion and 1/N C expansion, h q -> f q ml = 0.0025 GeV 3 , which is the 
value used in our scheme. As discussed in Ref. [14], the full expression (27) yields 
h q = (0.0015 ± 0.004) GeV 3 , i.e. a 200% uncertainty, if the errors of / 9jS and are treated 
as uncorrelated. The large error is due to a cancellation between the two terms in (27). 
As the parameter we need is actually h q /(2m q ), with m q not very well constrained (yet) 
from lattice calculations 3 and the correlation of the errors of f q>s and is not known, 
we feel that a total 250% uncertainty of h q / (2m q ) is slightly exaggerated and an artifact 
of the numerical cancellation. Instead, we work to leading order in the chiral expansion 
and set h q /(2m q ) = f q B , with B = m 2 /(2m (? ) = -2(0|gg|0)// 2 [29]. (0|gg|0), the 
quark condensate, is the order parameter of chiral symmetry breaking and known from 
QCD sum rules to have the value (0|gg|0) = (-0.24 ± 0.01) 3 GeV 3 . From this, one finds 
Bq = (1.6 ± 0.2) GeV [29], which, together with the error on f q , implies a total 15% 
uncertainty for the normalisation of the twist-3 DAs. This is the standard treatment of 
these terms in the framework of light-cone sum rules. 

3 LCSRs for Gluonic Contributions 

The key idea of light-cone sum rules is to consider a correlation function of the weak 
current and a current with the quantum numbers of the B meson, sandwiched between 
the vacuum and an rj or rj' state. For large (negative) virtualities of these currents, the 
correlation function is, in coordinate-space, dominated by distances close to the light- 
cone and can be discussed in the framework of light-cone expansion. In contrast to 
the short- distance expansion employed by conventional QCD sum rules a la SVZ [37], 
where non-perturbative effects are encoded in vacuum expectation values of local operators 
with vacuum quantum numbers, the condensates, LCSRs rely on the factorisation of 
the underlying correlation function into genuinely non-perturbative and universal hadron 
DAs 0. The DAs are convoluted with process-dependent amplitudes T H , which are the 

3 A recent unquenched calculation yields to = (m u + m d )/2 = (3.54+Q35) MeV at the scale /x = 2 GeV 
[36]. 
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analogues of the Wilson coefficients in the short-distance expansion and can be calculated 
in perturbation theory. Schematically, one has 

correlation function ~ ^ <g> <f> n . (28) 

n 

The expansion is ordered in terms of contributions of increasing twist n. The light-cone 
expansion is matched to the description of the correlation function in terms of hadrons by 
analytic continuation into the physical regime and the application of a Borel transforma- 
tion, which introduces the Borel parameter M 2 and exponentially suppresses contributions 
from higher-mass states. In order to extract the contribution of the B meson, one de- 
scribes the contribution of other hadron states by a continuum model, which introduces a 
second model parameter, the continuum threshold s Q . The sum rule then yields the form 
factor in question, / + , multiplied by the coupling of the B meson to its interpolating field, 
i.e. the B meson's leptonic decay constant /#. 

LCSRs are available for the B — > ir, K form factor / + to 0(a s ) accuracy for the 
twist-2 and part of the twist-3 contributions and at tree-level for higher-twist (3 and 4) 
contributions [26]. 

We define the B — > P form factors as 

{ p {M B (P+q) ) = { (2P + - sm . (29) 

Note that we include a factor 1 / \[2 on the right-hand side. This is to ensure that in the 
limit of SU(3)f symmetry and no 77-77' mixing j\ — /+. 

In the semileptonic decay B — > rf^lvi the form factor f$ (P = 77, 77') enters proportional 
to the lepton mass mf and hence is irrelevant for light leptons (/ = e,/i), where only /f 
matters. The semileptonic decay can be used to determine the size of the CKM matrix 
element \V u b\ from the spectrum 

rJV H 2 IT/ I 2 

where X(x) = (m 2 B + m 2 p — x) 2 — 4m 2 B m 2 p . Alternatively, as we shall see, the ratio 
of branching ratios B(B — > v(iv)jB(B — > rjiu) can be used to constrain the gluonic 
Gegenbauer moment B\. 

Our starting point for calculating is the correlation function 

KM = * f ^e^)(P(p)|T[^6](x)^(0)|0) (31) 

= K(q 2 ,p 2 B )(2p + q) l , + ... 

where jb = m b ui^ 5 b is the interpolating field for the B meson and p 2 B = (p + q) 2 its 
virtuality. For 



m 



1-P%> 0(A QCD m fe ), m 2 b -q 2 > 0(A QCD m b ), (32) 
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the correlation function (31) is dominated by light-like distances and therefore accessible 
to an expansion around the light-cone. The above conditions can be understood by de- 
manding that the exponential factor in (31) vary only slowly. The light-cone expansion is 
performed by integrating out the transverse and "minus" degrees of freedom and leaving 
only the longitudinal momenta of the partons as relevant degrees of freedom. The inte- 
gration over transverse momenta is done up to a cutoff, // IR , all momenta below which 
are included in a the DAs (fi n . Larger transverse momenta are calculated in perturbation 
theory. The correlation function is hence decomposed, or factorised, into perturbative 
contributions T and nonperturbative contributions <fi, which both depend on the longitu- 
dinal parton momenta and the factorisation scale /irR. The schematic relation (28) can 
then be written in more explicit form, including only two-particle DAs, as 

lW,p|) =J2 f duT^(u,q 2 , P 2 B ^ m )cf> n .p(u^ m ). (33) 

n J 

As Il + itself is independent of the arbitrary scale /zir, the scale-dependence of and 
(fi n must cancel each other. If there is more than one contribution of a given twist, they 
will mix under a change of /ira, and it is only in the sum of all such contributions that the 
residual /xir dependence cancels. This is what happens with the two-quark and two-gluon 
contributions to B — > r)('\ Eq. (33) is called a "collinear" factorisation formula, as the 
momenta of the partons in P are collinear with the P's momentum. Its validity actually 
has to be verified, which is done precisely by checking that the /zni dependence cancels. In 
Ref. [26] it has been shown that the above formula holds to 0(a s ) accuracy for two-quark 
twist-2 and -3 contributions. 

In calculating the correlation function, we use relation (8) between I?/'-*) and the QF 
basis states |% s ), so that 

n^ = i=(n«cos0-n^sin0), n£ = i= (n* sin <fi + n* cos 0) . (34) 

As the correlation function involves the current «7 M 6, IT* vanishes to leading order in a s 
and at 0(a s ) is due only to gluonic Fock states of the meson. IT^, on the other hand, 
receives contributions from both quark and gluon states. The quark contributions have 
been calculated in Ref. [26] for B — > n, including 0(a s ) corrections to twist-2 and -3 
contributions, and to tree-level accuracy for twist-4 contributions. The corresponding 
expressions yield IT+, with the replacement /„■ — > f q . 

In order to obtain the singlet contribution to IT^, one needs to calculate the diagrams 
shown in Fig. 3. The projection of the gluon fields onto the DA ifi 9 2 . p can be read off 
Eq. (15). The explicit formula is given in the appendix. We check the result by verifying 
the cancellation of the /xpR-dependent terms as described above. The relevant term in the 
quark Gegenbauer moment a 2 is 

K ~ W q F(pl, <f)a 2 (l + ^ f In , (35) 
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Figure 3: Feynman diagrams of the gluonic contributions. The double line denotes the b quark, 
the photon- like lines the currents in the correlation function IP^\ The first diagram is divergent, 
the other two are convergent. 

where F(p 2 B ,q 2 ) is a function of p 2 B and q 2 . The logarithmic terms in the convolution of 
the gluonic diagrams of Fig. 3 with ip^.p read 

One can easily convince oneself by expressing f q via Eq. (11) in terms of /jj and 
respectively, and inserting (35) into (34), that the renormalisation-group equation (21) is 
fulfilled. 

The final LCSR for /£ then reads 

e -m%/M> m 2 fB f±(<fi = P dse s/M^ lmYl P( SA 2^ (37) 
V2 Jml ^ 

with the sum-rule specific parameters M 2 , the Borel parameter, and s , the continuum 
threshold. 

4 Results and Discussion 

Let us now give the results for the form factors. As usual, we replace fs in the sum rule 
(37) by its QCD sum rule to 0(a s ) accuracy; this reduces the dependence of the results 
on mj = (4.80 ± 0.05) GeV. In Fig. 4 we plot /+(0) and /+ (0), respectively, as functions 
of the Borel parameter M 2 . The continuum threshold is chosen as so — 34.2 GeV 2 , 
which corresponds to the optimum s f° r the sum rule for [26]. The factorisation 
scale /xir is chosen as intermediate between m b and an intrinsic hadronic scale 1 GeV; 
following our earlier papers, we choose /z 2 R = m 2 B — m 2 . The dependence of on M 2 
is small in the Borel-window M 2 > 6 GeV 2 . We estimate the uncertainty in M 2 as the 
variation of the form factor in the interval M 2 G [6, 14] GeV 2 . In Fig. 4, we also show the 
dependence of the form factors on s by varying it by ±0.7 GeV 2 ; also this dependence 
is rather small. The central values of the most relevant hadronic input parameters are 
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Figure 4: [Colour online] /+(0) (left) and (0) (right) as a function of the Borel parameter 
M 2 and various choices of input parameters. Solid curves: central values of input parameters 
and so = 34.2 GeV 2 . Long-dashed (blue) curves: so varied by ±0.7 GeV 2 . Short-dashed (green) 
curves: 02(1 GeV) varied by ±0.15. Dash-dotted (red) curves: £>f varied by ±10. 
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Figure 5: [Colour online] /+(0)//^(0) as a function of the Borel parameter M 2 and various 
choices of input parameters. Solid (blue) line: central values of input parameters, which corre- 
sponds to /+(0)//^(0) = tan<rj = 0.814. Dash-dotted (red) curves: i?f varied by ±10. Short- 



dashed (green) curves: 6% (lGeV) varied independently: = 0.1, a 



0.4 and a 2 



0.4, 



0.1. 




q 2 [GeV 2 ] q 2 [GeV 2 ] 

Figure 6: [Colour online] /+(</ 2 ) (left) and /? (q 2 ) (right) as a function of the momentum transfer 
q 2 and various choices of input parameters. Solid curves: central values of input parameters and 
M 2 = 10 GeV 2 , s = 34.2 GeV 2 . Long-dashed (blue) curves: s varied by ±0.7 GeV 2 and M 2 
by ±4 GeV 2 . Short-dashed (green) curves: 02(1 GeV) varied by ±0.15, f q /f n by ±0.02 and 4> 
by ±1°. Dash-dotted (red) curves: -Bf varied by ±10. 
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mi, — 4.8 GeV, a 2 ,v (1 GeV) = 0.25 and B 9 = 0. As expected, /+(0) is not very sensitive 
to the singlet contribution parameter B 2 (red/dashed-dotted curves), but rather sensitive 
to the Gegenbauer moment a,2 (green/short-dashed curves). For /+(0), on the other 
hand, the dependence on B 2 is more pronounced than that of <i2- Varying all relevant 
parameters within their respective ranges, i.e. Am b = ±0.05 GeV, Aa 2 (lGeV) = ±0.15 
and AB 2 = ±20, as well as all twist-3 and twist-4 parameters within the ranges given in 
Ref. [30], we find 

fl(0) = 0.229 ± 0.005(M 2 ) ± 0.006(s ) ± 0.016(4') ± 0.007(Sf) ± 0.005(/„ 0) 

±0.011(T3) ±0.001(T4) ±0.007(/ B ,m 6 ) 
= 0.229 ± 0.024(param.) ± 0.011(syst.) , (38) 
/£(0) = 0.188 ± 0.004(M 2 ) ± 0.005(s ) ± 0.013(4'') ± 0.043(S|) ± 0.005(/„ 0) 

±0.009(T3) ± 0.005(T4) ± 0.006(/ B , m h ) 

= 0.188 ± 0.0025f ± 0.019(param.) ± 0.009(syst.) . (39) 

The entry labelled T4 also contains an estimate of the possible impact of the local twist-4 
two-gluon matrix elements in (19). For this estimate, we exploit the fact that the asymp- 
totic DA of the non-local generalisation of (19) is the same as for the twist-2 two-quark 
DA: 6u(l — u). 4 We then assume that the corresponding correlation function is the same 
as that for the leading conformal wave in the two-quark twist-2 contribution, i.e. the coef- 
ficient in the Gegenbauer moment ao = 1, and replace ao by (0\a s GG / (4n)\r] q ^) / \f q , s mf) . 
The factor 1 /m\ comes from the fact that this is a twist-4 effect and hence suppressed by 
two powers of mb with respect to the twist-2 contribution. This is only a rough estimate, 
of course, as the true spectral density will be different. The result in (39) shows that for 
small B 9 pa 2 both twist-2 and -4 two-gluon effects can indeed be of similar size. In this 
case, however, the total flavour singlet contribution to f_f_ will also be small, ~ 0.008. In 
the third lines, we have added all uncertainties from the input parameters (param.) in 
quadrature and the sum-rule specific uncertainties from M 2 and s (syst.) linearly. For 
/+ (0), we have displayed the dependence on B 9 2 separately. Our new result for /+(0) is, 
within errors, in agreement with our previous one, /+(0) = 0.275 ± 0.036, obtained in 
Ref. [26]. That for /+ (0) is new. Our results agree well with those obtained in Ref. [23], 
from perturbative QCD factorisation, /+(0) = 0.208 and f\ (0) = 0.171, including a 
rescaling by a factor \pl to bring their definition of the form factors into agreement with 
ours. We confirm the finding of Ref. [23] that the range of the singlet contribution to the 
form factor estimated in Ref. [14] is likely to be too large, unless B 2 assumes extreme 
values ~ 40. 

4 This follows from the general formula for asymptotic DAs, u 2jl_1 (l — u) 2 - 72-1 , with j = 1/2(1 + s) 
the lowest conformal spin of the operator, and I its canonical dimension, s the Lorentz-spin projection. 
For Gj_j_, one has I = 2 and s = [30]. 
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In Fig. 5 we plot the ratio f \ (0)//+(0) as a function of the Borel parameter. In the 
ratio, many uncertainties cancel, in particular that on As we have chosen B 2 = 
as central value, /+(0)//+(0) = tan0 = 0.814 exactly, see Eq. (34). The figure also 
illustrates the change of the result upon inclusion of a non-zero B 2 (red/dashed-dotted 
curves) . The ratio is actually rather sensitive to that parameter. While the dependence on 
a 2 largely cancels when a 2 and a 2 are set equal, there is a considerable residual dependence 
on a 2 — a 2 ^ (green/short-dashed curves). While \a 2 — \ = 0.3 as illustrated by 
these curves is rather unlikely, and would signal very large OZI-breaking contributions 
(recall that a 2 ^ a 2 or, equivalently, a\ ^ a 2 signals the presence of "wrong-flavour" 
contributions to the i] qs DAs and is set to in the QF mixing scheme), one should 
nonetheless keep in mind that moderate corrections of this type are not excluded and 
compete with the OZI-allowed corrections in B 2 . 

Let us now turn to the dependence of the form factors on q 2 . In Fig. 6 we show 
this dependence in the range < q 2 < 16 GeV 2 accessible by LCSRs. Again we dis- 
play in blue (by long-dashed curves) the dependence of /+ (q 2 ) on the sum-rule specific 
parameters M 2 and s , the green (short-dashed) curves illustrate the dependence on a 2 
and other parameters and the red (dash-dotted) ones that on B 2 . We give two different 
parametrisations of the form factors, in terms of a sum of two poles, the so-called BZ 
parametrisation as given in Ref. [26], and in terms of the BGL parametrisation based 
on analyticity of /+ in q 2 [38]. Both parametrisations are fitted to the LCSR results 
in the range < q 2 < 16 GeV 2 , and can then be used to extrapolate these results to 
<?max = ( m B — Tn v (.>)) 2 ', this is possible as both parametrisations include the essential 
feature of the B*{1) pole at q 2 = m 2 Bt ., m B * = 5.33 GeV, which governs the large-g 2 
behaviour of b — > u vector-current transitions close to q^^- 

The BZ parametrisation reads 

U(q 2 ) = /+(0) (- — Ijs- + — 2/ ?w7^ 27-W ) ' ( 4 °) 

\l-q 2 /m 2 B , (1 - q 2 /m 2 B „) (1 - aq 2 /m B ) J 

with the two shape parameters a, r and the normalisation /+(0). The BGL parametrisa- 
tion, on the other hand, is given by 

1 x 

«&^^ 9 » )f ' (41) 

{q 2 + - q 2 Yl 2 - {q 2 + - g 2 } 1 / 2 
{q% - Q 2 Y /2 + {q% - q 2 } 1/2 ' 

(g+ - q 2 )(V~q+ -q 2 - + Vq+- q 2 f l2 Wq+ -q 2 + Vq+- go) 
(V%+ VqT^Tiql-qlY" " ' 

(m B ±m v(l) ) 2 . (42) 

The "Blaschke" factor P(q 2 ) = z(q 2 , m B *) accounts for the B* pole, q 2 is a free parameter 
that can be chosen to attain the tightest possible bounds, and it defines z^q^ql) = 0. 



U(q 2 ) = 
with z(q 2 ,ql) = 

<f>(q 2 ,q 2 ) = 

and q\ = 
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Figure 7: /+ ' '(g 2 ) for central values of input parameters, fitted to the BGL parametrisation 
(41), for < k max < 9. 



One has \z\ < 1 for q 2 < (mg + m v (,)) 2 . In the following we choose g 2 , such that z(0, q 2 ,) = 



% 2 ) 



i.e. g 2 



14.14GeV for r] and 10.85 GeV for 77'. With these values, 



-z(q 2 _ 

becomes minimal: \z\ < 0.13 for rj and \z\ < 0.08 for 77'. The series in (41) provides a 
systematic expansion in the small parameter z, which for practical purposes has to be 
truncated at order k max . In this paper, we choose /c max = 3. 

The advantage of the BZ parametrisation is that it is both intuitive and simple: it 
can be obtained from the dispersion relation for / + , 



/fV) = — ^ 



U{<f 



+ 



m. 



7T 



(m B +m (,)) 



dt 
2 t 



lmff(t) 



ie 



(43) 



by replacing the second term on the right-hand side by an effective pole. However, it can- 
not easily be extended to include more parameters. The strength of the BGL parametri- 
sation, on the other hand, is that the dominant behaviour in q 2 close to the pole at m 2 B » is 
factored out and the remaining g 2 -dependence is organised as a Taylor-series in the small 
g 2 -dependent parameter z; the truncation of the series can be adjusted to the accuracy 
of the available input parameters. In Fig. 7 we plot j\ (q 2 ) parametrised a la BGL for 
< A: m ax < 9. Obviously, the parametrisations converge rapidly with increasing /c max and 
only differ at very large q 2 . The impact of this difference on the predicted branching ratio 
(30) is however only minor, as this region is phase-space suppressed. In the following, we 
choose k max = 3, which ensures that the total predicted branching ratio agrees within 1% 
with that obtained for k max = 9. 

In Tab. 1 we give the best-fit parameters for /+ in the BZ parametrisation, with the 
small effects of non-zero B 9 2 expanded linearly in that parameter. Tab. 2 contains the 
corresponding parameters for the BGL parametrisation with /c max = 3. Finally, in Fig. 8 
we show the dependence of the ratio of branching ratios R m i = B(B — > r]'ev)/B{B — > r\ev) 
on £>f . The advantage of this observable is that all hadronic effects are encoded in the 
form factors and that \V u b\ cancels. The blue (solid) curve corresponds to the branching 
ratios obtained from the central values of input parameters; the dependence of these 
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/+(0) 


a 


r 


V 
rf 


n 9 oi +0.018 

U.ZOl_ 020 

n 1 8Q+ - 015 _L U9 /+0.002N 

u.ioy_ i6 ~r r>2 \-0.w2) 


n sc; 1+0.183 

U.OOl_ 492 

fl +0-185 , R 9 /-0.006\ 
U.OOl_ 4g7 -|- r> 2 ^+0.008,/ 


n 41 1 -0.030 
u.^ii +0 119 

u -* li +0.122 ^ - D 2 1,-0.006,1 



Table 1: Parameters for the BZ parametrisation (40). The uncertainty contains all sources of 
error added in quadrature, except for rf, where the uncertainty in B 2 is approximated by a linear 
term. The upper (lower) terms represent the maximum (minimum) value of the form factor. 





V 


rf 


a 


0.0031 ±0.0003 


0.0018 ± 0.0002 ±0.00002£f 


ai 


-0.0090 =F 0.0034 


-0.0058 =F 0.0016 q= 0.00015f 


a 2 


0.0243 ±0.0172 


0.0174 ±0.0166T0.00015f 


a 3 


-0.0908 T- 0.0039 


-0.1189 =F 0.0218 ± 0.0016fif 



Table 2: Like Tab. 1, but for the BGL parametrisation (41) with /c max = 3. 

predictions on the cut-off in k is very small: the long-dashed (blue) curves illustrate the 
dependence on k max — 3 ± 1. On the other hand, also depends on a 2 7^ a\ . This 
dependence is shown by the red (short-dashed) curves. The conclusion is that large values 
of B%, \B 9 2 1 > 5, can be distinguished from the OZI-breaking parameter \a 2 —a 2 |, once an 
accurate experimental value of R m > is available, but that for smallish B 2 and unknown 
\a 2 — a 2 I only mutual constraints on these parameters can be extracted from the data. 
In this mentioned before, also twist-4 gluonic DAs can become important. 

To summarise, we have calculated the form factors of B — ■> rf'} semileptonic transitions 
from QCD sum rules on the light cone, including the gluonic singlet contributions. We 
have found that, as expected, these contributions are more relevant for /+ than for fl 
and can amount up to 20% in the former, depending on the only poorly constrained 
leading Gegenbauer moment B 9 2 of the gluonic twist- 2 distribution amplitude of rf'\ We 
also found that the form factors are sensitive to the values of the twist-2 two-quark 
Gegenbauer moments a 2 ,v which, given the uncertainty of independent determinations, 
we have set equal to a 2 . The ratio of branching ratios B(B — > r]'ev)/B{B — > r\ev) is 
sensitive to both a 2 and B 2 and may be used to constrain these parameters, once it is 
measured with sufficient accuracy. The extraction of \V u b\ from these semileptonic decays, 
in particular B — > rjeu, with negligible singlet contribution, although possible in principle, 
at the moment is obscured by the lack of knowledge of a 2 . We would also like to stress 
that, in the framework of the quark- flavour mixing scheme for the rj-rj' system as used 
in this paper, B — > rf'} transitions probe only the r] q component of these particles. The 
i] s component could be probed directly for instance in the b — > s penguin transition 
B s — > vf'}i + i~ , although such a measurement would also be sensitive to new physics in 
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Figure 8: [Colour online] The ratio of branching ratios R m ' = B{B — ► r}'ev)/B{B — > rjev) as a 
function of the singlet-parameter S| • Solid (blue) curve: central values of input parameters and 
BGL parametrisation with k max = 3; long-dashed (blue) curves: BGL parametrisations with 
&max varied by ±1. Short-dashed (red) curves: theoretical uncertainty of R m < for £?f = 0, for 
a™ (1 GeV) varied independently, as in Fig. 5. 



the penguin diagrams. 
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A Spectral Density of the two-gluon Contribution to 

/+ 

The contribution of the twist-2 two-gluon distribution amplitude to the correlation func- 
tions IT}, and iPj., Eq. (34), is given by 



2 

m\ S ~ Pb 



with 



pf (s) = Bla s f(m h - A= {59mt + 21g 6 - 63g 4 s - 19gV + 2s 3 

36V3 [s - q 2 f 

+m 2 s(164g 2 + 13s) - m 4 (82g 2 + 95s)} 
+ Bia s f[m b ^= iml ~f}* 2 ~ ml) {5m 4 + g 4 + 3g 2 s + s 2 - 5m 2 (g 2 + s)} 

x ( 2 ln^i-ln4V (A.l) 
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